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1. Introduction 

More than 40 years ago Ginibre [1] proposed three types of general Gaussian non- 
^ ' Hermitian matrix ensembles, those with complex, quaternion real and real entries. While 
he already gave the joint probability density (jpd) of eigenvalues and correlations for 
\0 ' the first two cases, the last case was harder to solve. It started with the presentation of 
^ ' the jpd by Lehmann and the author 25 years later [2]. Edelman rederived the jpd and 
O ■ determind the density of complex eigenvalues in 1997 [3]. In more recent publications 
Kanzieper and Akemann presented the complex correlations as Pfaffians [1], Sinclair 
^ ' derived a generating Pfaffian functional [5] and finally Forrester and Nagao were able 

■ to determine the real and complex correlations as Pfaffians with the help of orthogonal 
polynomials [6]. The last step should be considered as the solution of the problem of 
correlations of the real Ginibre ensemble. Here we give an alternative simple derivation 
and the generalisation to arbitrary real, complex or crossed correlations. 
In this paper we present the general n-point densities, all entries real or complex, as 
quaternion determinant governed by the simple kernel 

* 71=0 
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We restrict ourselves to even dimension A^. This is the only place where the matrix 
dimension occurs. The power expansion of ([1]) generates the inverse of an A^- 
dimensional submatrix {1 ^ k,l ^ N) of the infinite dimensional matrix 

Aki = j d^z, j d^Z2 :F{z,,z^)z\~^zi-^ (2) 

where T{zi, Z2) is a skew-symmetric complex measure reflecting the symplectic structure 
of the real Ginibre ensemble. A simple argument without use of orthogonal polynomials 
yields directly equation ([1]). Furthermore we discuss properties of l-and-2-point 
functions in more detail. 

The real asymmetric Gaussian ensemble has many applications in physics and social 
sciences as biological webs [Tj, neural networks [8], directed quantum chaos [9], financial 
markets [10]. The probability density of matrices Jij (1 ^ i,j ^ A^) we are talking 
about in this paper is simply oc exp(— Jfj/2). And the eigenvalues Zi of the matrix 
Jij are the zeros of the characteristic polynomial det{Jij — z) = and therefore real or 
pairwise complex conjugate. 



2. Symplectic Structure 

We start with the joint probability density of eigenvalues for the real Ginibre ensemble 
derived by Lehmann and Sommers [2] and calculate from it a generating functional for 
correlation functions [5l [6] : 

Z[f] = j d^^i ■ ■ ■ d^^N P{zi, 2:2, • • • , ZN)f{zi) . . . /(zn) . (3) 

P contains [y] + 1 pieces with R real eigenvalues {0 ^ R ^ N) and Q = {N — R)/2 
complex conjugate pairs in a chosen order. We extend here for convenience the 
eigenvalues to complex values, Zk = Xk + iyu- Thus for real eigenvalues, z = Xk, P{. ■ .) is 
concentrated on the real axis: i/k = 0. Afterwards we may symmetrize P{zi, Z2, ■ ■ ■ , zn) 
with respect to all permutations of Zk and integrate all variables Zk over the whole 
complex plane dividing the result by A^!. P{. . .) is proportional to the Vandermonde 
determinant 

l[{z, - z,) = det {z1-\ . . . , z'^-') (4) 

i>j 

which has a definite sign for the chosen order. Using the method of alternating variables 
[12] we may then integrate out ([2D and obtain a Pfaffian. The simplest way to derive 
it, is using Grassmannians, see the discussion at the end of this paper. For even A^ the 
result is 



Z[f] = U Pfaff (AfcO = M\Jdet{A>,i) (5) 
with a skew-symmetric matrix 

Aki = j d^z,d^Z2f{zr)f{z2)Hz,,Z2)z\-h'2^ (6) 
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where J-'{zi,Z2) is a skew-symmetric function [Ti^Zx^Z'^ = —J-'{z2, Zi)) of two complex 
variables zi, Z2 

T{zx.Z2) = e-(^?+^^')/2 [m\z, - Z2){Q{y,) erfc(yiv^) 

- 6(1/2) erfc(y2V2)} + 6{y,)6{y2){e{x2 - x^) - Q{x, - X2))] . (7) 

A/" is a normalization constant, which can be restored afterwards. It is related to the 
volume of the real orthogonal group [2J. Actually the variables zi, Z2 are considered here 
as two dimensional real vectors. This symplectic measure reveals two contributions, one 
from the real axis with |/i = 1/2 = and one from the complex plane off the real axis 
with S^{zi — Z2) = S{xi — X2)6{yi + j/2). The step functions 0(. . .) reflect the chosen 
order of eigenvalues. Functional derivatives of Z[f] with respect to f{z) at / = 1 give 
immediately all n-point densities: 



Riizi) 



(8) 



/.I Sf{z,)6f{z2) 

Equations ([3]) - ([7]) are completely equivalent to the joined density of eigenvalues 
obtained by Lehmann and Sommers [2] and rederived by Edelman [3]. We will use them 
to calculate the correlation functions (or ra-point densities, which should be distinguished 
from connected correlations). 
From ([5]) we obtain 

mf] 

with the kernel 



Z[f] j d^Z2 f{z2)H^l, Z2)K{Z2, z,) (9) 



1...N 



K{z2,z^) = J2Am4-'4^'- (10) 



kl 



At / = 1 we have A^i = A^i, K{z2, zi) = K{z2, zi) and Z[l] = 1. Thus 

Ri{zi) = [ (fz2Hzi.Z2)K{z2,zi) . (11) 



Ri{zi) contains two independent parts, one on the real axis and one in the complex 
plane. Comparing the latter with Edelman's expression for the density in the complex 
plane (not on the real axis) |3j we obtain (since z^ and the complex conjugate z^ can be 
considered as independent variables) 

, I...N 

K(z,_. .0 = 5: = ^ ^-.4-14-. . (12) 



^ n=0 kl 



Thus, remarkably, the inverse A^-dimensional skew-symmetric matrix A^^ has a very 
simple tridiagonal form and the bulk part and the real-axis part of the complex density 
Ri{z) are intimately related via the kernel K{z2^ zi). Equation (|T2l) can independently 
be checked by calculating A^^ directly from ^ and ([7]). Now, not only we recover 
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Edelman's bulk part of Ri{z) [3], but also we obtain the density on the real axis first 
obtained by Edelman, Kostlan and ShubfTT] 

Ri{z) = R^\z) + 6{y)R^{x) (13) 



with 



R?iz)=e-^^'-y'^^-^ [ due--'\y\J^\^ (14) 



\y\V2 



and 



-oo 



^ / \ /I N—2 / /\ ^ 



2v^ ^ n! 



Writing 



^-'^ I i\n °? N-2 

-XX' ^ ' [xx'r I , 



n=0 



(16) 



we see that expression (fTS!) is positive for even ^ 2, for which it was derived. For 
odd there is a correction term, which ensures positivity. Equation (fT6|) makes it easy 
to derive large A^-expansions by saddle-point integration. The simplest results for large 
are the Girko circle law 



i?f(z) ~ -e(yiv- 1^1) (17) 

and a constant density on the real axis [6] 

/?f (x) ~ -^e(v^ - kl) . (18) 



This is in contrast to the semicircle density for a Gaussian Hermitian matrix (Wigner 
semicircle law). A new type of behaviour is expected near the edge. The average number 
of real eigenvalues Nfi{N) is in agreement with Edelman [3]: 



1 



...X V2 f dt t^/\l - 1^-') bN , , 

^n(^)-^ + -J (i-t)3/.(i+t) -Vv 







For calculating the next correlation functions one has to take functional derivatives of 
([9]), for which we use 

SK{z2,Zi) = J d^Z^5f{z-i) j d'^Z4,f{z4)J^{z3,Zi) ■ 

K{z2,z,)K{zs,z,)-K{z2,zs)K{z,,z,)'^ . (20) 

This equation is enough to generate all higher functional derivatives. Equation (!20l) 
implies for the 2-point density 

R2{Zl,Z2) = J^iZi,Z2)K{z2,Zi) + j d"^ Z3 j d"^ Z^J^iZi, Zs)J^{z2, Z4) ■ 

{K{zs, z,)K{z^, Z2) + K{z3, Zi)K{z2, z,) - K{zs, Z2)K{z,, Zi)} . (21) 
At the end of this paper we will write this as a quaternion determinant. 
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3. 2-point densities 

Let us single out in R2{zi, Z2) the nonsingular bulk part 2:2), that contains no 

(5-contributions and describes correlations in the complex plane off the real axis ^ 

R2{zi,Z2) = 4 sgn(?/i) erfc(||/i|\/2) sgn(t/2) erfc(|?/2| v^) ■ 

-K{z^, Z2)K{Z2, ^1) + K{z^, Z2)K{z-2, Zl)} (22) 

For small ^1,-22 this behaves as 16\yiy2\{x1 + Xg) > 0. For (-21,^1) — » (^2,^2) the 
correlation i?^(zi,2;2) behaves oc \zi — 2:2^. On the other hand, the most singular 
part of R2{zi,Z2) on the real axis is 6{yi)6{y2)R2 {xi, X2) where i?^(xi,X2) describes 
correlations of real eigenvalues and is given by: 

R^{xi,X2) = e-(^?+^2)/2 I sgn{x2 - Xi)K{x2, Xi) + 

dx3 J dxi e~*''^3+^4)/2 ggj2(^3;g _ nf.^'^ sgn(a;4 — X2) ■ 
■ [K{x'i,xi)K{xi,X2) + K{x'i,Xi)K{x2.,xi) - K{x'i,X2)K{xi,xi)\} . (23) 

Both expressions fl22|) and fl23l) are contained in results of [6j. We proceed now to derive 
scaling forms of the correlations in the large- iV limit. Since the densites ( ITTl) and ( fTSll 
in the complex plane and on the real axis are constant, the scaled correlation functions 
are simply obtained by taking the N 00 limit of equation (fT2i) 

K(z2,^i)^f^e^^^^. (24) 



Obviously the symplectic measure function jF(zi,z2) (JZ!) does not depend on A^. This 
implies for the 1-point densities 



R^{z) ~ y|e2^'|i/| erfc(|y|v^) 



(25) 



i?f (x) ^ . (26) 



and 



Both are independent of x and only for \y\ ^ 1, i.e. far from the real axis, Ri{z) goes 
to ( |T71) . This is generally valid: far from the real axis the correlations go to Ginibre's 
(i.e. the Ginibre ensemble for general complex matrices [1]). For example the bulk part 
i?^(zi, Z2) of the 2-point correlation function for large is given by 



R'i{zi,Z2) ~ R^Azi)R'{{z2)-^Uym 



4yil/2 

+ ((xi - X2f + (2/1 - y^f)e-^^--^-^'-^y^+y^)' 

-((xi - X2f + (2/1 + y^f)e-^^^-^^^'-^y--y^r^ (27) 

Far from the real axis only the Ginibre result 

i?2^(zi,Z2)^^(l-e-l^-^^l') (28) 
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survives. Finally we compute the asymptotic form of the real 2-point correlations 
R^{x,, X,) ^ ^ f 1 - e-(— )0 + ^^±^e-^ erfc ( . (29) 



27r V / V V2 

This is different from Wigner-Dyson correlations, but shows {3=1 level repulsion. For 
small distances we have i?f^(xi,a;2) — for large distances we have Rf{xi,X2) — 



27r \^ 



^ 2V2-K 

1 ~ ^(xl-x2)'^ j ) ^ much faster decay of connected correlations than for GOE (the 
Gaussian orthogonal ensemble), more similar to those for the Ginibre ensemble in the 
complex plane. 



4. n-point densities 

Finally we observe that the correlation (120!) can be written as a quaternion determinant 
of a selfdual matrix (related to a Pfaffian of a skew-symmetric matrix) and can be 
generalized for n-point densities {k,l = 1,2, ... ,n) 



j (Pz T{zk,z)K{z,zi) K{zk,zi) 

-T{zk, Zl) J <Pz K{zk, z)T{z, Zl] 



(30) 



RniZi, ...,Zn) = Qdet 

with 

J^{zi, Z2) = T{zx, Z2) - J (fz(fz'J^{zi, z)K{z, z)T{z , z^) . (31) 

The sign of the quaternion determinant is given by the diagonal elements 
Rx{z\), . . . ,Rn{zn) (see ffTTj) ). The term R\{z}) ■ R^iz^) ■ ■ ■ Rniz^^ appearing in the 
expansion of the quaternion determinant just gives the asymptotics of Rn{zi, . . . ,Zn) 
for large separation. 

For an alternative derivation one may write the generating functional Z[l + u] as a 
Gaussian integral over Grassmannians ( 

Z[l+u] = J DCexp(^-^ J d'z J d'z'az)az'mz,z')il + uiz))il + uiz'))^ (32) 

with Grassmannian fields ^{z) with {^{z)^{z'))u=o = K{z, z!) and expand Z[l +m] using 
the fermionic Wick theorem to obtain all correlation functions. The entries in (l30l) 
are the building elements of the diagrammatic expansion. A similar expansion for the 
quaternion determinants proves the claim fl30l) . One crucial point is that although the 
quaternion determinant is a square root, it is an analytic function of the entries. 



5. Conclusion 



In conclusion, without use of orthogonal polynomials we have given a transparent deriva- 
tion of the correlations for the real Ginibre ensemble, which are governed by one simple 
symplectic kernel K[^z, z') describing a deep connection between real and complex cor- 
relations. Real eigenvalues, which meat at a point on the real axis, continue to move 
along the imaginary direction repelling each other without changing the real coefficients 
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of their characteristic polynomial. Calculations have been done for even A^. The algebra 
for odd is more complicated. However one may conjecture, that the n-point densities 
for odd N can be obtained by analytic continuation from even preserving positivity. 
Moreover we have discussed some properties of l-and-2-point functions in detail. A lot 
of further properties can be discussed in future. 

After completing most of this paper, I have seen some recent papers [H O [6], where 
many of the results reported here are partly contained. However this paper contains in 
addition some generalisations and tries to make transparent the beauty of the structure 
of the real Ginibre ensemble, which also may be useful in many fields of sciences. 
The author acknowledges support by SFB/TR12 of the Deutsche Forschungsgemein- 
schaft. 
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